Plectoneme formation in twisted fluctuating rods

The mechanics of DNA supercoiling is a subject of crucial importance to uncover the mechanism and kinetics of several enzymes. It is therefore being investigated using several biochemical and biophysical methods including single molecule experimental techniques. An interesting problem within this realm is that of torsional buckling and plectoneme formation in DNA as it is simultaneously put under tensile and torsional stress. Analytical solutions to this problem are difficult to find since it involves nonlinear kinematics and thermal fluctuations. In this paper we use ideas from the Kirchhoff theory of filaments to find semi-analytical solutions for the average shape of the fluctuating DNA under the assumption that there is no self-contact. The basic step in our method consists of combining a helical solution of the rod with a non-planar localizing solution in such a way that the force, moment, position and slope remain continuous everywhere along the rod. Our solutions allow us to predict the extension vs. linking number behavior of long pieces of DNA for various values of the tension and temperature. An interesting outcome of our calculations is the prediction of a sudden change in extension at buckling which does not seem to have been emphasized in earlier theoretical models or experiments. Our predictions are amenable to falsification by recently developed single molecule techniques which can simultaneously track the force–extension as well as the torque–rotation behavior of DNA

Equilibrium and Kinetics of DNA Overstretching Modeled with a Quartic Energy Landscape

AbstractIt is well known that the dsDNA molecule undergoes a phase transition from B-DNA into an overstretched state at high forces. For some time, the structure of the overstretched state remained unknown and highly debated, but recent advances in experimental techniques have presented evidence of more than one possible phase (or even a mixed phase) depending on ionic conditions, temperature, and basepair sequence. Here, we present a theoretical model to study the overstretching transition with the possibility that the overstretched state is a mixture of two phases: a structure with portions of inner strand separation (melted or M-DNA), and an extended phase that retains the basepair structure (S-DNA). We model the double-stranded DNA as a chain composed of n segments of length l, where the transition is studied by means of a Landau quartic potential with statistical fluctuations. The length l is a measure of cooperativity of the transition and is key to characterizing the overstretched phase. By analyzing the different values of l corresponding to a wide spectrum of experiments, we find that for a range of temperatures and ionic conditions, the overstretched form is likely to be a mix of M-DNA and S-DNA. For a transition close to a pure S-DNA state, where the change in extension is close to 1.7 times the original B-DNA length, we find l ≈ 25 basepairs regardless of temperature and ionic concentration. Our model is fully analytical, yet it accurately reproduces the force-extension curves, as well as the transient kinetic behavior, seen in DNA overstretching experiments

A continuum model for the growth of dendritic actin networks

Polymerization of dendritic actin networks underlies important mechanical processes in cell biology such as the protrusion of lamellipodia, propulsion of growth cones in dendrites of neurons, intracellular transport of organelles and pathogens, among others. The forces required for these mechanical functions have been deduced from mechano-chemical models of actin polymerization; most models are focused on single growing filaments, and only a few address polymerization of filament networks through simulations. Here we propose a continuum model of surface growth and filament nucleation to describe polymerization of dendritic actin networks. The model describes growth and elasticity in terms of macroscopic stresses, strains and filament density rather than focusing on individual filaments. The microscopic processes underlying polymerization are subsumed into kinetic laws characterizing the change of filament density and the propagation of growing surfaces. This continuum model can predict the evolution of actin networks in disparate experiments. A key conclusion of the analysis is that existing laws relating force to polymerization speed of single filaments cannot predict the response of growing networks. Therefore a new kinetic law, consistent with the dissipation inequality, is proposed to capture the evolution of dendritic actin networks under different loading conditions. This model may be extended to other settings involving a more complex interplay between mechanical stresses and polymerization kinetics, such as the growth of networks of microtubules, collagen filaments, intermediate filaments and carbon nanotubes

Entropically Driven Motion of Polymers in Nonuniform Nanochannels

In nanofluidic devices, nonuniform confinement induces an entropic force that automatically drives biopolymers toward less-confined regions to gain entropy. To understand this phenomenon, we first analyze the diffusion of an entropy-driven particle system. The derived Fokker-Planck equation reveals an effective driving force as the negative gradient of the free energy. The derivation also shows that both the diffusion constant and drag coefficient are location dependent on an arbitrary free-energy landscape. As an application, DNA motion and deformation in nonuniform channels are investigated. Typical solutions reveal large gradients of stress on the polymer where the channel width changes rapidly. Migration of DNA in several nonuniform channels is discussed

Particle diffusion in active fluids is non-monotonic in size

We experimentally investigate the effect of particle size on the motion of passive polystyrene spheres in suspensions of Escherichia coli. Using particles covering a range of sizes from 0.6 to 39 microns, we probe particle dynamics at both short and long time scales. In all cases, the particles exhibit super-diffusive ballistic behavior at short times before eventually transitioning to diffusive behavior. Surprisingly, we find a regime in which larger particles can diffuse faster than smaller particles: the particle long-time effective diffusivity exhibits a peak in particle size, which is a deviation from classical thermal diffusion. We also find that the active contribution to particle diffusion is controlled by a dimensionless parameter, the Peclet number. A minimal model qualitatively explains the existence of the effective diffusivity peak and its dependence on bacterial concentration. Our results have broad implications on characterizing active fluids using concepts drawn from classical thermodynamics.Comment: 5 Figure

Effect of supercoiling on formation of protein-mediated DNA loops

DNA loop formation is one of several mechanisms used by organisms to regulate genes. The free energy of forming a loop is an important factor in determining whether the associated gene is switched on or off. In this paper we use an elastic rod model of DNA to determine the free energy of forming short 50–100 basepair, protein mediated DNA loops. Superhelical stress in the DNA of living cells is a critical factor determining the energetics of loop formation, and we explicitly account for it in our calculations. The repressor protein itself is regarded as a rigid coupler; its geometry enters the problem through the boundary conditions it applies on the DNA. We show that a theory with these ingredients is sufficient to explain certain features observed in modulation of in vivo gene activity as a function of the distance between operator sites for the lac repressor. We also use our theory to make quantitative predictions for the dependence of looping on superhelical stress, which may be testable both in vivo and in single-molecule experiments such as the tethered particle assay and the magnetic bead assay

Force steps during viral DNA packaging?

Biophysicists and structural biologists increasingly acknowledge the role played by the mechanical properties of macromolecules as a critical element in many biological processes. This change has been brought about, in part, by the advent of single molecule biophysics techniques that have made it possible to exert piconewton forces on key macromolecules and observe their deformations at nanometer length scales, as well as to observe the mechanical action of macromolecules such as molecular motors. This has opened up immense possibilities for a new generation of mechanical investigations that will respond to such measurements in an attempt to develop a coherent theory for the mechanical behavior of macromolecules under conditions where thermal and chemical effects are on an equal footing with deterministic forces. This paper presents an application of the principles of mechanics to the problem of DNA packaging, one of the key events in the life cycle of bacterial viruses with special reference to the nature of the internal forces that are built up during the DNA packaging process

A statistical mechanics framework for constructing non-equilibrium thermodynamic models

Far-from-equilibrium phenomena are critical to all natural and engineered systems, and essential to biological processes responsible for life. For over a century and a half, since Carnot, Clausius, Maxwell, Boltzmann, and Gibbs, among many others, laid the foundation for our understanding of equilibrium processes, scientists and engineers have dreamed of an analogous treatment of non-equilibrium systems. But despite tremendous efforts, a universal theory of non-equilibrium behavior akin to equilibrium statistical mechanics and thermodynamics has evaded description. Several methodologies have proved their ability to accurately describe complex non-equilibrium systems at the macroscopic scale, but their accuracy and predictive capacity is predicated on either phenomenological kinetic equations fit to microscopic data, or on running concurrent simulations at the particle level. Instead, we provide a framework for deriving stand-alone macroscopic thermodynamics models directly from microscopic physics without fitting in overdamped Langevin systems. The only necessary ingredient is a functional form for a parameterized, approximate density of states, in analogy to the assumption of a uniform density of states in the equilibrium microcanonical ensemble. We highlight this framework's effectiveness by deriving analytical approximations for evolving mechanical and thermodynamic quantities in a model of coiled-coil proteins and double stranded DNA, thus producing, to the authors' knowledge, the first derivation of the governing equations for a phase propagating system under general loading conditions without appeal to phenomenology. The generality of our treatment allows for application to any system described by Langevin dynamics with arbitrary interaction energies and external driving, including colloidal macromolecules, hydrogels, and biopolymers